Astronomy & Astrophysics manuscript no. SS02rv2 


February 1, 2008 


(DOI: will be inserted by hand later) 





Spin-down of Relativistic Stars with Pliase Transitions and PSR 

J0537-6910 

N. K. Spyrou and N. Stergioulas 



m 
o 
o 

(N 
D 
l> 

> 
O 
00 

m 
o 

O 



Oh- 
I ■ 

o 



X 



Astronomy Department, Aristoteleion University of Thessaloniki, 541.24 Thessaloniki, Macedonia, Greece 
email: spyrou@astm.auth.gr, niksterg@astro.auth.gr 

Received / Accepted 

Abstract. An abstract siiould be given 



Using a highly accurate numerical code, we study the spin 
down of rotating relativistic stars, undergoing a quark decon- 
finement phase transition. Such phase transitions have been 
suggested to yield an observable signal in the braking index 
of spinning-down pulsars, which is based on a "backbending" 
behaviour of the moment of inertia. We focus on a particular 
equation of state that has been used before to study this be- 
haviour, and find that for the population of normal pulsars the 
moment of inertia does not exhibit a backbending behaviour. 
In contrast, for supramassive millisecond pulsars a very strong 
backbending behaviour is found. Essentially, once a quark core 
appears in a spinning-down supramassive milUsecond pulsar, 
the star spins up and continues to do so until it reaches the insta- 
bility to collapse. This strong spin-up behaviour makes it eas- 
ier to distinguish a phase transition in such pulsars: a negative 
first time-derivative of the rotational period, P < 0, suffices and 
one does not have to measure the braking index. In the spin- 
up era, the usually adopted spin-down power law fails to de- 
scribe the evolution of the angular velocity. We adopt a general- 
relativistic spin-down power law and derive the equations that 
describe the angular velocity and braking index evolution in 
rapidly rotating pulsars. We apply our numerical results to the 
fast young pulsar J0537-6910 in SNR N157B, which has been 
suggested to have (if spun down by magnetic dipole radiation 
only) an extremely small initial spin period. The inclusion of 
a quark-hadron phase transition can yield a significantly larger 
initial spin period of 6ms (in our example), which is in bet- 
ter agreement with theoretical expectations. Finally, we suggest 
that the frequent rate of glitches in PSR J0537-6910 could be 
related to the fact that it is the fastest Crab-Uke pulsar, so that a 
pure quark core may have formed recently in its lifetime. 

1. Introduction 

The properties of the highest-density region of matter in the in- 
terior of relativistic stars have been studied to date by several 
different theoretical models, which yield very different macro- 
scopic properties of compact stars. Many new interactions have 



been proposed to occur at the highest densities attained, rang- 
ing from the occurrence of hyperons, to quark deconfinement, 
to the formation of kaon condensates and H-matter (see Weber, 
2001 ; Heiselberg 2002, for recent reviews). Observational con- 
straints on pulsar properties are currently still too weak to al- 
low to distinguish which (if any) of all the proposed theoreti- 
cal descriptions of the interior of relativistic stars is the correct 
one. A solution to the problem of determining the correct equa- 
tion of state in the central region of a compact star must come 
from new, more accurate observations of their properties and 
any new observational method that will help in this direction is 
more than welcome. 

One new observational method, that could be interesting 
in the case that the equation of state (EOS) features a decon- 
finement phase transition to quark matter, has been proposed 
by Glendenning, Pei and Weber (1997), (hereafter GPW). 
Specifically, when a rapidly rotating pulsar spins down, its cen- 
tral density increases with time. At a certain central density, 
a mixed quark-hadron phase can appear (or can already be 
present in the pulsar). For some chosen parameters of the equa- 
tion of state and of the baryonic mass of the pulsar, a pure quark 
core appears when the pulsar reaches very high densities. In 
GPW. a particular spin-down sequence was studied, that in the 
nonrotating limit reached nearly at the maximum mass allowed 
by the chosen EOS. For this sequence, it was noticed that, when 
the pure quark core appears, the star undergoes an brief era of 
spin-up (brief, compared to the pulsar's lifetime) and its brak- 
ing index shows an anomalous behaviour (it becomes singu- 
lar at a certain rotation rate). This change in the braking index 
has been proposed in GPW to be potentially observable. The 
anomalous behaviour of the spin evolution is traced back to a 
"backbending" of the moment of inertia (as a function of angu- 
lar velocity). 

Other authors (Heiselberg & Hjorth-Jensen 1998, 
Chubarian et al. 2000) have further investigated the ob- 
servational consequences in the models proposed in GPW and 
in Glendenning and Weber (2000). In all these studies the 
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relativistic slow-rotation approximation is used. In the present 
work, we re-investigate the deconfinement phase-transition 
of spinning-down pulsars, using fully relativistic, rapidly 
rotating models and find that, when the quark-core appears (in 
a sequence of models similar to the one used in GPW), the 
backbending in the moment of inertia is generically absent for 
normal pulsars, but very strong for supramassive millisecond 
pulsars. Furthermore, our numerical results suggest that in the 
limiting region between normal and suprammasive pulsars, the 
backbending behaviour is sensitive to truncation errors of the 
applied numerical scheme. The moment of inertia is a very 
sensitive function of angular velocity and must be computed 
with high accuracy. Truncation errors, such as a limited num- 
ber of points in the EOS table, or the use of the slow-rotation 
approximation can be large enough to systematically affect the 
computation of the moment of inertia. 

The case of supramassive millisecond pulsars is very in- 
teresting, as the spin-up era lasts from the onset of the quark- 
core appearence to their collapse to black holes, i.e. for their 
remaining lifetime. For supramassive millisecond pulsars cre- 
ated via accretion-induced spin-up, the spin-up era could last 
for the most part of their lifetime, if a quark-hadron phase tran- 
sition occurs. Thus, the first time-derivative of the rotational 
period can be used as a tool for observing quark-hadron phase 
transitions and one does not have to rely on the measurement 
of the braking index. We emphasize that the extended spin-up 
era in supramassive millisecond pulsars is not restricted for the 
case of a quark-hadron phase transition, but it can also occur 
for other types of phase transitions, such as a phase transition 
to a pion condensate. An example is EOS M (Pandharipande 
and Smith, 1975) for which an extended spin-up with angu- 
lar momentum loss (for supramassive, constant baryonic mass 
sequences) has been found in Cook, Shapiro and Teukolsky 
(1994). On the other hand, the absence of spin-up in the cur- 
rently known millisecond pulsar population favours the ab- 
sence of large phase transitions or implies that the EOS is such 
that millisecond pulsars created via accretion-induced spin-up 
are probably not supramassive. 

In previous studies of spinning-down pulsars that were lim- 
ited to the slow-rotation approximation (including up to 0{Q?) 
rotational effects) expressions for the evolution of angular ve- 
locity and braking index were derived (see Glendenning, 1997). 
However, we find that these expressions are incomplete, as 
changes in the gravitational mass of the star (other than changes 
in the kinetic energy) were not taken into account in the as- 
sumed spin-down law. Here, we adopt a general-relativistic ver- 
sion of the spin-down law, which dictates that the energy lost 
in the form of e.g. magnetic dipole radiation or gravitational 
waves is not lost to the expense of the kinetic energy only, but 
to the expense of the total mass-energy of the star. Hence, mak- 
ing use of a relation due to Bardeen (1970), which expresses 
the first law of thermodynamics along sequences of constant 



After this work was completed, we learned of a new study using 
rapidly rotating relativistic stars, with a numerical code related to ours, 
(Cheng, Yuan & Zhang, 2002), in which the backbending behaviour 
of the moment of inertia is shown to be sensitive to the assumed prop- 
erties of the crust. 



baryonic mass, we derive expressions for the angular velocity 
and braking index evolution, which are valid for any rotation 
rate. In addition, our adopted spin-down law is also valid dur- 
ing a spin-up era, while it's Newtonian, slow-rotation limiting 
expression (which is usually adopted in the literature) is not 
valid during spin-up. 

As an application of our numerical investigations, we com- 
pute a possible value for the initial period of pulsar J0537- 
6910, assuming that the high-density EOS exhibits a quark- 
deconfinement phase transition and that the star spins down 
only due to magnetic dipole radiation. Without the assump- 
tion of a phase transition, previous studies suggested that the 
initial spin of the pulsar would be extremely small (for a brak- 
ing index of n = 3) (Marshall et al. 1998). Here we show that 
the presence of the phase transition can yield a more reason- 
able initial period for this pulsar, without having to assume a 
braking index smaller than n = 3. In our example, we obtain 
a significantly larger initial period of 6ms, which is in better 
agreement with theoretical expectations. Furthermore, this pul- 
sar has displayed a large number of giant glitches (Gotthelf et 
al. 2001). We propose that these could be due to the appear- 
ance of a pure quark core in the center of the star, which results 
in an increased rate of contraction of the star during its spin- 
down. In contrast, rapidly rotating pulsars do not exhibit such 
giant glitches, when their central density is less than the criti- 
cal density for the appearance of the quark core. The fact that 
pulsar J0537-6910 displays such glitches more frequently than 
other Crab-like^ pulsars may be related to the fact that it is the 
most rapidly spinning of all Crab-like pulsars and thus a pure 
quark core could have appeared more recently in its lifetime. 

The outline of the paper is as follows: In Sect. |2l we de- 
scribe the chosen equation of state and its numerical refine- 
ment, while in Sec.|3]we obtain the equilibrium properties of 
a specific sequence of spinning-down compact stars. The be- 
haviour of the moment of inertia and gravitational mass along 
the sequence is described in Sec.l^and thermodynamic consis- 
tency is used as an independent test in Sec.|5l In Sec.|6lexpres- 
sions for the evolution of the angular velocity and the brak- 
ing index are derived, and in Sec. the initial period of PSR 
J0537-6910 is estimated. Finally, in Sec.|8j we comment on the 
expected event rate of the considered signal. We conclude with 
a discussion of our results in Sec.|9l 

2. Equation of State and Numerical Refinement 

The EOS used in GPW features a mixed quark-hadron phase, in 
between pure hadron and pure quark phases. The mixed phase 
is obtained by applying the Gibbs criteria for two-component 
systems, abandoning local (but not global) charge-neutrality 
(Glendenning 1992). This yields a second-order phase transi- 
tion between each region. However, it has been suggested (see 
Heiselberg & Hjorth- Jensen, 1999) that the mixed phase can 
exist only if the interface tension between quarks and hadrons 
is not too large. But, on the other hand, the interface tension 
may not be a free parameter, but it could have the freedom to 

^ In this paper, "Crab-like" only refers to the young age of a pulsar 
and to its observational association with a supernova remnant. 



N. K. Spyrou and N. Stergioulas: Spin-down of Relativistic Stars with Pliase Transitions 



3 



adjust (via a minimization prinicple) to the nature of the two 
phases (Glendenning 2002), see Christiansen, Glendenning & 
Schaffner-BieHch (2000) for a relevant example. 

The specific EOS we use is the quark deconfinement EOS 
given in Table 9.2 of Glendenning (1997) (same as in GPW) in 
which, however, we replace the pure quark phase by an an- 
alytic expression (see discussion below). The pure confined 
hadronic phase exists below baryon number density 0.246fm^^ 
while pure deconfined quark matter exists above 0.862fm^^ 
(lfm=10^'^cm). In between these two phases, a mixed quark- 
hadron phase exists (see Fig.[n which displays the pressure P 
vs. energy density e relation for this EOS). 
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Fig. 1. Tabulated points of pressure P vs. energy density e for 
the quark deconfinement EOS. 



When computing nonrotating models with the original tab- 
ulated EOS, using the 2-D code by Stergioulas & Friedman 
(1995), we noticed that for central energy densities near the 
central energy density where the pure quark phase appears, 
the gravitational mass vs. central energy density plot shows 
an unusual behaviour (see Fig. |2j- There is a local maximum 
followed by a local minimum and the mass increases again, 
before it reaches the maximum mass allowed by the general- 
relativistic radial instability (see Shapiro & Teukolsky, 1983). 
Such a behaviour is not present, for exactly the same tabu- 
lated EOS, when using a simple 1-D code that integrates the 
well-known radial TOV equations of structure for non-rotating 
models. Thus, the oscillatory behaviour of the results obtained 
with the 2-D code for nonrotating models is a numerical ef- 
fect. In the 2-D code, the angular behaviour of all variables 
is expanded as a sum of Legendre polynomials Pi (/j) of argu- 
ment jj = cos 9. At discontinuities, such as the surface of the 
star (where e.g. the velocity profile is discontinuous), the ac- 
curacy of the code is reduced due to the Gibbs phenomenon 
(see Nozawa et al. 1998). Similarly, when the pure quark phase 
appears at the center of the star, discontinuities appear in the 
numerical solution, as the EOS features an abrupt change of 
slope (see Fig.[0. In addition, as the EOS table only contains 
a limited number of points, the interpolation methods (such as 
4-point Lagrange interpolation or Hermite polynomial interpo- 
lation) that are used in our code, do not avoid an oscillatory 



behaviour in the interpolated quantities, near the discontinuity 
in the slope of the EOS. It appears that, in the 2-D code, such 
an interpolation error enhances the Gibbs phenomenon at the 
center of the star, when the pure quark core appears, resulting 
in the oscillatory behaviour in Fig.|2l 

To solve this numerical problem, we replace the part of the 
EOS that represents the pure quark phase by an analytic ex- 
pression. In fact, we find that this part can be represented rather 
accurately by the following linear relation: 

P[dynes/cm^] = 1.85240 X 10^^ (1) 
+2.885209 X \Q^^{t/c^ - 1.69100 x lO'^) 

where P is pressure, e is energy density and c is the speed of 
light (in CGS units). In this way, the tabulated EOS is used only 
in the region below the pure quark phase and no interpolation is 
done through the phase transition, thus eliminating the problem 
of the oscillatory behaviour For nonrotating models, this re- 
fined tabulated/analytic EOS now reproduces the same mass vs. 
central energy density relation with the 2-D code, as the orig- 
inal tabulated EOS with the 1-D code (see Fig.|2j- Notice that 
for this comparison we avoid interpolating through the region 
where the pure quark appears, by choosing points either before 
or after the phase transition only. This increases the accuracy 
of computed models (as compared to interpolating through all 
points in the table) even with the 1-D code. 

Thus, we have verified that replacing the pure quark region 
of the original tabulated EOS by the analytic expression (|2j 
eliminates numerical problems in the 2-D code, while still rep- 
resenting the same EOS (with sufficiently high accuracy). 




1.413 
1.6e+15 



1.8e+15 2e+15 

zjc (g/cm^) 



2.2e+15 



Fig. 2. Gravitational Mass M vs. central energy density for 
nonrotating models, computed by various methods: using the 
original tabulated EOS in the 2-D code produces an oscilla- 
tory behaviour (dotted line). In contrast, using the refined tab- 
ulated/analytic EOS in the 2-D code (solid line) gives results 
in agreement with the original tabulated EOS in the 1-D code 
(when interpolation through the phase transition is avoided) 
(filled circles). 
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3. Spin-down Sequence 

We study several evolutionary sequences of pulsars, that differ 
from each other in the assumed constant baryonic mass. For 
a baryonic mass of Mq = 1.555Mp) we obtain a sequence that 
terminates at the maximum mass nonrotating model allowed 
by the EOS. This sequence is the limiting sequence between 
normal pulsars and supramassive pulsars and is very similar 
to the sequence studied in GPW. We assume that a compact 
star with a baryonic mass of Mq is rapidly rotating. The com- 
pact star subsequently spins down due to e.g. magnetic braking 
and/or gravitational radiation emission, and the sequence ends 
either at the nonrotating limit (for normal pulsars), or the pul- 
sar encounters the axisymmetric instability limit to black hole 
collapse (see e.g. Cook, Shapiro and Teukolsky, 1984, for ex- 
amples of such sequences for various realistic EOSs). The se- 
quences are constructed using a highly accurate numerical code 
(rns, Stergioulas & Friedman, 1995) that computes exact nu- 
merical solutions of rapidly rotating stars in general relativity, 
without any approximation (apart from the truncation errors of 
the numerical scheme). Solutions are obtained on a 2-D numer- 
ical grid, with typical grid-size of 800x400 or 1200x400 (radial 
X angular) points. Each model is defined by its central energy 
density and polar to equatorial axis ratio. Apart from the re- 
quirement of high accuracy for each individual model, essential 
to the construction of a constant-baryon number sequence, is 
the accuracy with which a model of a given baryonic mass is lo- 
cated. To this end, we construct various models by keeping the 
central energy density fixed, varying only the axes ratio and use 
a robust root-finding method to locate a model of given bary- 
onic mass within a specified accuracy. The root-finding method 
of our choice is Ridders' method (Ridders, 1979), since the pro- 
cedure is guaranteed to stay within the initial brackets (avoid- 
ing unphysical parameters) and since it has a relatively high 
order of convergence (\/2)- In practice, about 10-12 individual 
models are required for the root-finding method to converge to 
a model of given baryonic mass, with a relative accuracy of 
10^^ (the accuracy in computing the baryonic mass for an in- 
dividual model is of the order of 10^^, but this does not limit 
the root-finding method to converge to a specified value with 
higher relative accuracy). 



4. Moment of Inertia and Gravitational Mass 

We compute several physical parameters of the equilibrium 
models, the moment of inertia /, gravitational mass M, angular 
momentum J and angular velocity Q. being the most relevant 
for our present study. We find that the behaviour of the mo- 
ment of inertia depends on the baryonic mass of the sequence. 
Along a normal sequence with Mq = I.SSIMq, which is very 
close to the maximum mass nonrotating limit, but still below it, 
the change in the moment of inertia, as the star spins down, is 
shown in Fig.|3] which plots / as a function of Q. (solid Une). 
Using our refined tabulated/analytic EOS, described above, the 
behaviour of the moment of inertia is seen to be regular, i.e. / is 
a single-valued function of Q.. We notice that an apparent back- 
bending appears with our 2-D code for this sequence, when the 
EOS in the original tabulated form (and interpolation through 



all points in the table) is used (dashed line in Fig. |3}- This is 
an indication that the backbending of the moment of inertia is 
sensitive to the available points in the EOS table and/or to trun- 
cation of rotational effects on the structure of the relativistic 
star. 




3000 



Fig. 3. Moment of Inertia / vs. angular velocity £2, along 
the normal spin-down sequence with Mq = I.SSIM©. The 
original tabulated EOS produces an apparent backbending in 
the moment of inertia (dashed line) while the refined tabu- 
lated/analytic EOS yields the evolution of the moment of in- 
ertia without a backbending behaviour. 



Another important equilibrium property that we monitor 
along the equilibrium sequence is the gravitational mass. In 
Fig. |4] the gravitational mass M is plotted as a function of 
central energy density 8c . The solid line corresponds to the 
refined tabulated/analytic EOS and shows a monotonous de- 
crease of the gravitational mass during spin-down. This normal 
behaviour corresponds to the monotonous behaviour of the mo- 
ment of inertia in Fig.|3l In contrast, when one uses the original 
tabulated form of the EOS, the gravitational mass appears to 
oscillate as a function of angular velocity. This oscillatory be- 
haviour is not correct, as it is the result of numerical truncation 
errors, due to the interpolation through the phase transition in 
the original tabulated EOS. 

The oscillatory behaviour obtained with the original tabu- 
lated EOS deserves some further attention. As we will show, 
if it were the correct behaviour, it would correspond to a part 
of a sequence being unstable to quasi-radial perturbations: In 
this sequence, before the quark core appears, the equilibrium 
models are stable, while the star loses kinetic energy T and an- 
gular momentum J, and the radius decreases only due to the 
decrease in rotational flattening of the star during spin down. 
When the pure quark core appears at the center of the star, 
the radius starts decreasing significantly, as the EOS becomes 
softer Looking at models with larger central densities along the 
equilibrium sequence, the angular velocity f2, angular momen- 
tum J, kinetic energy T , absolute value of gravitational bind- 
ing energy |W| and gravitational mass M, all increase. Thus, 
a star cannot continue evolving along this uniformly rotating 
sequence to the expense of the star's available mass-energy. 
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Fig. 4. Gravitational Mass M vs. central energy density Ec along 
the normal spin-down sequence with Mq = I.SSIM©. Using 
the original tabulated EOS produces two turning points (that 
would correspond to an unstable region) (dashed line) while 
the refined tabulated/analytic EOS produces a monotonous be- 
haviour. 

In Fig. 13 (dashed lined) a turning point appears and the grav- 
itational mass increases with increasing central density. This 
is a typical case of the onset of the instability to the funda- 
mental axisymmetric, quasi-radial mode of oscillation, in ro- 
tating relativistic stars, as shown by Friedman, Ipser & Sorkin 
(1988) (hereafter FIS). The turning-point criterion for the on- 
set of instability was first proved in FIS for a constant-angular- 
momentum sequence. Along such a sequence, an extremum in 
the M vs. Ec curve separates the regions of stable and unstable 
stars, i.e. 



)j 



<^ onset of instability, 



(2) 



for a sequence of uniformly rotating stars obeying a one- 
parameter P = P(e) equation of state. Such a criterion was then 
proved by Cook, Shapiro & Teukolsky (1992) to also exist for 
sequences of constant baryonic mass Mq: 



<^ onset of instability. 



(3) 



Mo 



The instability to the fundamental quasi-radial mode of os- 
cillation is first a secular one, i.e. it proceeds on a viscous 
timescale. It has been shown, however, that a dynamical in- 
stability sets in near the secular instability point (see Shibata, 
Baumgarte & Shapiro 2000). Thus, if the star would evolve 
along the dashed curve in Fig. |3 reaching the secular insta- 
bility point, it would first evolve away from the uniformly- 
rotating equilibrium sequence, on a viscous timescale, driven 
by magnetic braking torques and/or gravitational radiation re- 
action and could then also become dynamically unstable, in 
which case it would suddenly contract. Eventually, the star 
would bounce at higher densities, when reaching again a stable 
branch in the equilibrium sequence (after the second turning 
point in Fig.l^. 

We find that, when computing sequences of constant bary- 
onic mass, the behaviour of the gravitational mass can be used 



Fig. 5. Moment of Inertia / vs. angular velocity Q., along sev- 
eral spin-down sequences. Normal sequences (Mq < l.SSSM©) 
do not show a backbending behaviour, while for supramassive 
sequences, such as the Mo = 1 .6M0 sequence, the backbending 
(spin-up) behaviour is dominant. 
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Fig. 6. Close-up of Figure [S] A small backbending behaviour 
first occurs for the limiting sequence of Mq = 1 .SSSM© that 
divides normal from supramassive stars and terminates at the 
maximum-mass nonrotating model allowed by the given EOS. 



as a guide to distinguish whether the sequence has been com- 
puted accurately: a behaviour such as the apparent instability 
described above, indicates a strong influence of numerical er- 
rors in the computation. 

Looking at other constant baryonic mass sequences, with 
baryonic mass between Mq = I.SSIM© and Mq = l.SSSM© 
(see Fig. |5j, we see that the moment of inertia gradually ap- 
proaches a vertical behaviour (dl/dQ. as the limiting 
sequence that divides normal from supramassive stars is ap- 
proached. A very small backbending (spin-up) behaviour is 
first observed for the limiting sequence of Mq ~ 1 .555M0 (and 
clearly not for e.g. Mq = 1.554M©, see Fig. |6}. For larger 
baryonic masses (supramassive sequences) the backbending 
behaviour becomes stronger and for e.g. a baryonic mass of 
Mq = I.6M0 it becomes the dominant behaviour, i.e. once the 
quark core appears, the star only spins-up with angular momen- 
tum loss, for the rest of its lifetime. 
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To summarize, for the particular EOS that we have studied 
in this paper, normal pulsars do not feature a spin-up era during 
phase transition, while supramassive pulsars only spin up after 
the phase transition. The limiting sequence that divides normal 
from supramassive pulsars shows a very small era of spin-up 
and appears to be the dividing line between stars that can spin 
up and stars that only spin down. 

5. Thermodynamic Consistency 

In the previous sections we demonstrated that the refined tab- 
ulated/analytic EOS produces a constant baryonic mass se- 
quence with equilibrium properties that are in sharp contrast 
to the equilibrium properties computed on the basis of the orig- 
inal tabulated EOS. In this section we will use an independent 
check to verify that the equilibrium properties obtained with the 
refined EOS are indeed the physically acceptable ones, while 
the equilibrium sequence obtained with the tabulated EOS suf- 
fers from numerical errors. 

The independent check of the accuracy of the computed 
equilibrium properties is provided by the use of a relation that 
was first proved by Ostriker and Gunn (1969) in Newtonian 
theory and then derived in general relativity by Bardeen (1970). 

Along a sequence of uniformly rotating models of constant 
baryonic mass, changes in M and J are related by 



c-'dM = Q.dJ, 



(4) 



which can be regarded as an expression of the first law of ther- 
modynamics for such sequences. 
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Fig. 7. Error indicator X, constructed with the use of the ther- 
modynamic consistency relation (see text) for the original and 
for the refined EOSs, along the sequence of rotating mod- 
els undergoing the quark-hadron phase transition, with Mq = 
1.551M0. The refined EOS gives significantly more accurate 
models in the region where the pure quark core appears. 

We have evaluated the above relation numerically along the 
obtained sequences of equilibrium models (using sufficiently 
many individual models to ensure an accurate evaluation of the 
numerical derivative). We define an error indicator, X, as 



1 



Q.dJ 



(5) 



The evaluation of X along the evolutionary sequence of Mq — 
1 .55 IM0 (Fig.0, shows that the above relation is satisfied with 
an accuracy of ~ 1 % for models with central energy densities 
smaller than the energy density at which the pure quark-core 
appears. After the appearance of the pure quark core, numerical 
errors result in huge values of X, when using the original tabu- 
lated EOS. With the use of the refined tabulated/analytic EOS, 
these errors are reduced significantly to an acceptable level and 
diminish as the resolution is increased. 

The above check confirms that the refined EOS indeed pro- 
duces physically acceptable results, while the original tabulated 
EOS produces results that are dominated by numerical error in 
the region in question. 



6. Behaviour of the Observable Braking Index 

If a compact star is detected as a pulsar and one can measure 
the first two time-derivatives Cl, Q. of its angular velocity, then 
one can define an observational braking index 



--n(D.) = 



(6) 



In the Newtonian, slow-rotation Umit, the spin-down of a pulsar 
is usually modeled in the form of a power law (see Shapiro & 
Teukolsky. 1983), namely assuming that only kinetic energy is 
lost and that the rate of loss is proportional to some power of 
the angular velocity of the star, i.e. 



dT 
dt 



(7) 



where K < and a > are usually assumed to be real con- 
stants. With these assumptions, the braking index is equal to 
n = a — 1 . For example, it is assumed that, for magnetic brak- 
ing, t ^ n^, which yields an expected braking index n — 2>, 
while, for gravitational wave emission, t ^ which yields 
« = 5. For slowly rotating pulsars, the above considerations 
are, in most cases, appropriate. However, for rapidly rotating 
pulsars one has to take into account the rotational flattening of 
the star To a first approximation, this can be done consider- 
ing rotational effects up to order (9(f2^). Glendenning (1997) 
gives the rotationally corrected Ci and braking index n(i2) for a 
sequence of uniformly rotating stars, assuming the spin-down 
law of Equation 0. As we will show here, the expressions 
given in Glendenning (1997) are incomplete, in the sense that 
the derivation is not fully consistent to but misses addi- 

tional contributions of the same order. 

The energy lost in the form of electromagnetic or grav- 
itation radiation is not only to the expense of the star's ki- 
netic energy (which would be the case only in the 0{£l) slow- 
rotation approximation) but to the expense of the star's total 
mass-energy (gravitational mass). In the Newtonian limit, this 
principle has been applied to the spin-down of rapidly rotat- 
ing neutron stars by Finn & Shapiro (1990) and by Spyrou & 
Stergioulas (2001). In relativistic stars, the total mass-energy is 



Mc^ = Mqc^ + U - 



-T + W, 



(8) 



N. K. Spyrou and N. Stergioulas: Spin-down of Relativistic Stars with Pliase Transitions 



7 



(see Friedman, Ipser & Parker 1986) where, additionally, U is 
the internal energy. We thus adopt the following natural gener- 
alization of Equation in general relativity: 



dt 



(9) 



The kinetic energy in relativity is defined as T = (1 /2)JQ., 
so that 



dT 
dt 



1 



,dM dO. 



2 \^ dt 



dt 



-J 



(10) 



Expressing this in terms of the moment of inertia, which is de- 
fined as in Friedman et al. (1986), / = and with the help 
of Equation (|4}, we find 



dT 
'dt 



,dM l^ndl 



dt 



dt 



(11) 



Equation shows that, in a rapidly rotating star, under the 
above assumptions, the time derivatives of the kinetic energy 
and gravitational mass differ by an 0{Q.^) term that is propor- 
tional to the time-derivative of the moment of inertia. Notice 
that this expression is exactly valid for any rotation rate (even 
for stars rotating at the mass-shedding limit) and not only in an 
0{Q?') approximation. The relations derived in Glendenning 
(1997) for the braking index and time-derivative of the angular 
velocity in the 0{Q?) approximation, are missing the contribu- 
tion from the 0{Q?) term on the r.h.s. of Equation il 1> . 



where /' = dl/dQ. and 



n(n) = a-l 



2I'Q. + I"Q?- 
I + I'Q. 



(13) 



We have verified that Equation (I13> is the correct one, by nu- 
merically computing the braking index as obtained above and 
comparing it to the direct computation of the braking index, 
using the definition Equation The two numerical results 
agree well even close to the mass-shedding limit. In contrast, 
the expression for n{D.) in Glendenning (1997) shows a large 
disagreement with the result obtained from the definition of 
n(n) for rapidly rotating stars. This difference between the two 
expressions is significant also in the case of supramassive se- 
quences. 

In Fig. Is] the observational braking index n{Q.) is shown 
as a function of angular velocity (computed along the nor- 
mal spin-down sequence with MO = l.SSlMp), assuming that 
the spin-down is only due to magnetic braking. The computed 
«(£2) is roughly n ~ 3 before the pure quark core appearance. 
During the phase transition, the braking index rises to roughly 
« ~ 10 and then decreases to a value of roughly n ~ —2, before 
returning to the canonical value of n ~ 3. This behaviour is 
qualitatively and quantitatively different from the singular be- 
haviour in the braking index obtained with our code when one 
uses the original tabulated EOS (see inset in Fig.|8}. 
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Fig. 8. Observational braking index n{Q.) along the spin-down 
sequence with Mo = I.SSIMq. The original tabulated EOS 
yields a singular behaviour at two nearby values of the angular 
velocity (see inset). The refined EOS does not yield a singular 
behaviour, but only a much weaker increase and corresponding 
decrease (with respect to the canonical value of n = 3) at the 
pure quark core appearance. 

With the above assumptions, it is easy to show that 
Equations (5.18) and (5 .20) in Glendenning ( 1 997) are replaced 
by 



/ 



(12) 



7. Initial Period of PSR J0537-6910 

PSR J0537-6910 was discovered in SNR N157B by Marshall 
et al. (1998). It has a spin period of only 16ms, the fastest of 
all Crab-like pulsars, and a period derivative of P — 5.13 x 
10^'^ss^'. Both the characteristic age z = P/{2P) and the age 
estimate for SNR N157B suggest that the pulsar is 5,000 years 
old. An upper limit on the age of the pulsar based on observa- 
tions of //a-emitting gas in N157B (Chu et al. 1992), is 20,000 
years. Marshall et al. show that, if the braking index is n = 3 
and the pulsar age is 5,000 years, then the pulsar must have 
been born with an extremely small (even sub-millisecond) pe- 
riod. Such small initial periods in pulsars associated with su- 
pernova remnants have not been detected to date. There are 
theoretical arguments, based e.g. on the possible occurrence of 
gravitational-wave induced instabilities in rapidly rotating rela- 
tivistic stars (see e.g. Stergioulas 1998, Andersson & Kokkotas 
2001 for reviews) that suggest a pulsar initial period (shortly 
after the compact star is born) of at least several ms. How can 
this theoretical expectation be compatible with the observed 
properties of PSR J0537-6910 ? One possible solution would 
be that the actual braking index for this pulsar is closer to n = 2 
(for the upper limit on age of 20,000 years, a very low braking 
index of « ~ 1 is required). 

Here we show that an initial period of several ms (6 ms 
in our example) can be obtained for a braking index of n = 
3, if the appearance of a pure quark core is considered. We 
use the same EOS and the MO = 1.551M0 constant baryonic 
mass sequence mentioned in previous sections and numerically 
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integrate Equation ( I12t . We find that the age of the pulsar can 
be expressed in the following integral form 



1 f^I + Q.dI/dQ. 
where £2, is the initial angular velocity, or, equivalently, 



1 /■■' 



age 



(14) 



(15) 



where 7, is the initial angular momentum. For numerical in- 
tegration, it is more convenient to introduce a dimensionless 
central energy density ic = (ec/c^)/(10'^g/cm^) as the inte- 
gration parameter, and evaluate 



'age 



1 



(16) 



The result of the numerical integration is shown in Fig.|9] 
which displays the age of the pulsar as a function of the central 
energy density. For a current spin of 16ms and a pulsar age of 
5,000 years, the initial spin is obtained to be 6ms. Thus, even 
an n = 3 braking index can yield an initial spin that is in agree- 
ment with theoretical expectations, provided the occurrence of 
a pure quark core is taken into account. In the above exam- 
ple, a particular choice of EOS and baryonic mass was made. 
Different choices would lead to different initial spin period esti- 
mates. However, the qualitative effect of the presence of a pure 
quark core is to increase (with respect to an EOS without a 
phase transition) the computed initial spin period. 
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Fig. 9. Computed age of PSR J0537-6910 as a function of its 
central energy density. For a current age of 5,000 years, the ini- 
tial spin period is 6ms, significantly larger than in the absence 
of a phase transition. 



8. Event Rate 

An important parameter for the successful detection of a phase- 
transition signal in the braking index of spinning-down pulsars 
is the event rate, which is determined by several important fac- 
tors. Let's assume that pulsars are born with an initial period 



of only a few milliseconds. Since normal pulsars have a mag- 
netic field strength of ~ 10 '^G, they would quickly spin down 
to much larger periods due to magnetic dipole radiation. Thus, 
if the baryonic mass of a newly-born pulsar is such that the 
phase transition appears at a short rotational period, then the 
phase transition signal will not last for an extended period of 
time to have a realistic chance of being detected. Notice that 
no young. Crab-like pulsar with period shorter than 16ms has 
been observed, which is indicative of the very small event rate 
that should be expected for such an observation. 

On the other hand, if pulsars are born with rotational pe- 
riods larger than, say, 6ms, then a phase transition can occur 
during their lifetime only if their baryonic mass falls within 
an extremely small range of values. To illustrate this point, we 
plot in Figure [^the constant baryonic mass sequence consid- 
ered in Sec. |3] and the corresponding sequence of nonrotating 
models. The vertical dashed line specifies the central density 
at which the pure quark core appears. It crosses the constant- 
baryonic-mass sequence at a model with rotational period of 
6ms and gravitational mass of \A\9Mq,. Notice that the gravi- 
tational mass of a nonrotating star with central energy density 
equal to that of the quark core appearance (1.691 x lO'^g/cm^) 
is 1.414M0. Thus, a pure quark core can appear in a spinning 
down pulsar with rotational period less than 6ms only if its 
mass falls within AM ^ O.OO25M0 of the value of 1.4165M0 
(the the sequence considered in this example). The probability 
for this to happen is obviously negligible. 

The above arguments show that the only observed popu- 
lation remaining, in which one could hope for a significant 
event rate for the considered signal, are old millisecond pul- 
sars, spun-up by accretion, with low magnetic field and spin- 
down rate. The actual event rate will then depend on the mass- 
distribution for this population and on the precise value of the 
central energy density at which the quark core appears. None 
of these parameters are currently known with sufficient accu- 
racy to compute a reliable value for the event rate. In fact, the 
observational determination of the braking index is itself a dif- 
ficult task and no braking index has been measured for any of 
the known millisecond pulsars, to date. However, should a mil- 
lisecond pulsar become supramassive, due to mass-accretion, 
then there is a realistic chance for detecting the presence of a 
quark core in observations of the first time-derivative, P, of the 
rotational period. Once the quark core appears and if the pulsar 
is sufficiently massive, then P <0 for the rest of the pulsar's 
lifetime. 



9. Discussion 

Observational constraints for the properties of compact stars 
are currently still too weak to determine the correct equa- 
tion of state of high-density matter The presence of a 
quark/hadron deconfinement phase transition has been sug- 
gested by Glendenning, Pei and Weber (1997) (GPW) to be 
observable in the braking index of pulsars. In this method, the 
observational signal comes from a spin-up era, originating from 
a backbending behaviour of the moment of inertia during the 
pure quark core appearance. 
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Fig. 10. Gravitational mass vs. central energy density Ec for 

the nonrotating sequence and for the spin-down sequence of 
constant baryonic mass Mq — I.SSIM©. The dashed vertical 
line indicates the central energy density at which a pure quark 
core appears. Slowly rotating normal pulsars can go through 
the phase transition, only if the gravitational mass falls within 
an extremely narrow range of values. 

We have investigated the deconfinement phase-transition of 
spinning-down pulsars, using fully relativistic, rapidly rotating 
models. We find that, for normal pulsars, the star does not go 
through a spin-up era and the pure quark core appears without 
a backbending behaviour in the moment of inertia. The back- 
bending behaviour first occurs for the Umiting sequence of con- 
stant baryon mass that terminates at the maximum mass nonro- 
tating model allowed by the EOS. In contrast, for sequences 
with larger mass (supramassive sequences) the backbending 
behaviour becomes dominant and pulsars only spin up, after 
a quark core appears. 

From our numerical results, it is evident that the moment 
of inertia along a spin-down sequence is very sensitive to the 
appearance of phase transitions and must be computed with 
great care. In particular, we emphasize that low resolution and 
the restriction to a limited number of points in the EOS table 
(with interpolation through all points) can easily spoil the cor- 
rect computation of the behaviour of the moment of inertia, at 
least in the present 2-dimensional code used in this paper. 

The expected event rate of a phase-transition signal in the 
pulsar braking index is negligibly small for the population of 
normal pulsars. For the population of old millisecond pulsars, 
a reliable estimate for the event rate cannot be obtained, as it 
involves several uncertain factors and, in addition, the obser- 
vational determination of braking indices in milUsecond pul- 
sars has not become technically possible, yet. However, if mil- 
Usecond pulsars become supramassive, then the presence of 
a quark core could be detected by observing a negative first 
time-derivative of the rotational period. On the other hand, 
the absence of any such observation in the current popula- 
tion of known millisecond pulsars indicates that either large 
phase transitions do not occur, or that the EOS is sufficiently 
stiff to prevent accreting relativistic stars from beconung supra- 
massive. A detailed study of phase transitions in supramas- 
sive models of relativistic stars, constructed with various EOSs, 



could provide interesting constraints on the equation of state at 
very high densities. 

In order to describe the spin down of a pulsar due to mag- 
netic dipole or gravitational radiation, we adopt a general- 
relativistic version of the usual power-law assumed in the slow- 
rotation limit. The energy loss is assumed to be to the expense 
of the total mass-energy of the star and not simply to the ex- 
pense of the star's kinetic energy. In this way, we derive new ex- 
pressions for the evolution of the angular velocity and braking 
index during spin down. As an application, assuming magnetic 
dipole braking only, the initial spin period of PSR J0537-6910 
is obtained to be 6ms (for the EOS and spin-down sequence 
considered in this paper). The presence of the deconfinement 
phase transition can thus give a more reasonable initial period 
for this pulsar, compared to previously obtained estimates. 

PSR J0537-6910 exhibits the highest rate of large glitch 
events (6 in 3 years) (see Gotthelf et al. 2001) of all known 
Crab-like pulsars. At the same time, it is the most rapidly rotat- 
ing Crab-hke pulsar. If a pure quark core appears at small spin 
periods, then the quark-core fraction in such pulsars continu- 
ously changes during spin down, most probably contributing 
to the occurrence of large glitches. The fact that PSR J0537- 
6910 exhibits the highest glitch rate could be related to the fact 
that the phase transition may have happened more recently in 
its lifetime (according to the spin-down sequence used here) 
and its spin-down behaviour is still influenced by the grad- 
ual change of the quark core fraction during spin down. This 
suggestion has to be examined in more detail, as it could be a 
generic feature of phase transitions in the high-density EOS. 
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